https://github.com/kzajac/Mownit2018/blob/master/Mownit_lab8.ipynb
using QuadGK using Polynomials using Plots using Roots